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[57] ABSTRACT 

A method for predicting the SEU susceptibility of a 
standard-cell D-latch using an alpha-particle sensitive 
SRAM, SPICE critical charge simulation results, and 
alpha-particle interaction physics. A technique utilizing 
test structures to quickly and inexpensively characterize 
the SEU sensitivity of standard cell latches intended for 
use in a space environment. This bench-level approach 
utilizes alpha particles to induce upsets in a low LET 
sensitive 4-k bit test SRAM. This SRAM consists of 
cells that employ an offset voltage to adjust their upset 
sensitivity and an enlarged sensitive drain junction to 
enhance the cell’s upset rate. 

3 Claims, 12 Drawing Sheets 
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METHOD FOR CHARACTERIZING THE UPSET 
RESPONSE OF CMOS CIRCUITS USING 
ALPHA-PARTICLE SENSITIVE TEST CIRCUITS 

5 

ORIGIN OF INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA Contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the Contractor has elected to retain title. 10 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

The present application is a continuation-in-part of 
U.S. patent application Ser. No. 07/672,705, filed on 15 
Mar. 19, 1991, now U.S. Pat. No. 5,331,164. 


TECHNICAL FIELD 


The present invention relates to bench-level assess- 
ment of digital IC sensitivity to SEU upsets and more 
specifically to a method for predicting the SEU suscep- 
tibility of a standard-cell D-latch using an alpha-particle 
sensitive SRAM, SPICE critical charge simulation re- 
sults, and alpha-particle interaction physics. 

BACKGROUND ART 


20 


25 


The trend in electronic components has lead to the 
development of integrated circuits (ICs) that have 
higher operating frequencies and require less power per 
gate. This has been achieved in part by a reduction in 30 
the feature sizes that make up the transistors and wires 
that comprise the ICs. As a consequence, the amount of 
charge that is required to store a bit of information is 
about one-tenth of a pico-coulomb (pC), or about 
600,000 electrons. Because the amount of charge is so 35 
small, cosmic rays can deposit enough charge in an IC 
to cause a single-event upset (SEU) and alter the state of 
the memory bit. 

The importance of cosmic rays on the performance of 
integrated circuits in a space environment is evident in 40 
the upset rates of various satellites and spacecraft. For 
example, the Tracking and Data Relay satellite experi- 
ences a single event upset per day, which must be cor- 
rected from the ground. Such adverse experiences have 
caused a re-design of spacecraft, such as the Galileo 45 
spacecraft. The characterization of digital ICs to heavy 
ion induced upsets is essential for qualifying their use in 
critical systems to be used on spacecraft. This charac- 
terization usually requires evaluations at an ion source, 
such as a cyclotron. Such evaluations are time consum- 50 
ing, expensive and error prone. In recent years, laser 
pulses have been proposed as substitutes for the heavy 
ion sources. However, the laser simulations are limited 
by metal layers that frequently block the laser pulses 
from SEU-sensitive nodes and by the complexity of the 55 
ion-photon calibrations. The solution to the single event 
upset problem continues to be important, as the com- 
plexity of spacecraft grows, the size of integrated cir- 
cuits decrease and as space systems are designed with 
circuits fabricated at non-radiation hardened foundries. 60 

Thus it would be highly desireable to have a method 
for bench-level characterization of an IC such as a 
CMOS standard-cell D-latch, without requiring the use 
of large and expensive heavy ion sources. 

SUMMARY OF THE INVENTION M 

The purpose of this invention is to provide a tech- 
nique utilizing test structures to quickly and inexpen- 


2 

sively characterize the SEU sensitivity of standard cell 
latches intended for use in a space environment. This 
bench-level approach utilizes alpha particles to induce 
upsets in a companion device, a low LET sensitive 4-k 
bit test SRAM. This SRAM consists of cells that em- 
ploy an offset voltage to adjust their upset sensitivity 
and an enlarged sensitive drain junction to enhance the 
cell's upset rate. 

The invention comprises a method for predicting the 
SEU susceptibility of a standard-cell D-latch using an 
alpha-particle sensitive SRAM, SPICE critical charge 
simulation results, and alpha-particle interaction phys- 
ics. Measurements were made on a 1.6-jitm n-well 
CMOS 4k-bit test SRAM irradiated with an Am-241 
alpha-particle source. A collection depth of 6.09 fim 
was determined using these results and TRIM computer 
code. Using this collection depth and SPICE derived 
critical charge results on the latch design, an LET 
threshold of 34 Mev cm 2 /mg was predicted. Heavy ion 
tests were then performed on the latch and an LET 
threshold of 41 MeV cm 2 /mg was determined. 

OBJECTS OF THE INVENTION 

It is therefore a principal object of the present inven- 
tion to provide a method for quickly and inexpensively 
characterizing the single-event upset (SEU) sensitivity 
of standard cell latches intended for use in a space envi- 
ronment. 

It is another object of the invention to provide a 
method for predicting the SEU susceptability of a 
standard-cell D-latch using an alpha-particle sensitive 
SRAM, critical charge simulation results and alpha-par- 
ticle interaction physics. 

It is still another object of the invention to provide a 
bench-level technique for characterizing CMOS stand- 
ard-cell D-latches using alpha-particle sensitive test 
circuits to obviate heavy ion sources which are expen- 
sive and time consuming to use. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The aforementioned objects and advantages of the 
present invention, as well as additional objects and ad- 
vantages thereof will be more fully understood herein- 
after as a result of a detailed description of a preferred 
embodiment when taken in conjunction with the fol- 
lowing drawings in which: 

FIG. 1 is a schematic drawing of a test SRAM cell; 

FIG. 2 is a schematic drawing of a D-latch in the zero 
state and illustrating reverse biased SEU-sensitive drain 
diodes; 

FIG. 3 is a layout drawing of a standard-cell latch; 

FIG. 4 is a graphical illustration of critical charge 
versus offset voltage for both an SRAM and a latch; 

FIG. 5 is a graphical illustration of upset rate versus 
offset voltage for an SRAM in response to a 4.6 /xCi 
Am-241 source; 

FIG. 6 is a graphical illustration of the relationship 
between measured upset rate, the spontaneous density 
function, the source spectra and the SRAM critical 
charge response; 

FIG. 7 is a graphical illustration of a test SRAM 
differential upset rate response versus offset voltage; 

FIG. 8 is a schematic illustration of the path of alpha- 
particles from an alpha-particle source through the 
SRAM; 

FIG. 9 is a graphical illustration of the D-latch alpha- 
particle response; 
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FIG. 10 is a graphical illustration of a heavy ion test 
performed on a latch for purposes of comparison with 
the inventive method herein; 

FIG. 11 is a graphical illustration of the SRAM criti- 
cal charge response for Node V2 of FIG. 1 calculated 5 
using SPICE-2; 

FIG. 12 is a graphical illustration of the number of 
flipped cells for alpha particles and protons shown rela- 
tive to the spontaneous bit flip response for a store time 
of one second; 10 

FIG. 13 is a schematic illustration of the paths of 
protons from a proton source to the SRAM; 

FIG. 14 is a graphical illustration of the SRAM upset 
probability distribution for alpha particles and protons 
shown relative to the spontaneous bit flip response; 15 
FIG. 15 is a graphical illustration of charge deposited 
by 0.55 MeV protons in silicon which is used to calcu- 
lated the over-layer thickness, 8X3; and 
FIG. 16 is a graphical illustration of charge deposited 
by 1.0 MeV protons in silicon which is used to calcu- ^ 
lated the collection depth, 8X4. 

DETAILED DESCRIPTION OF A PREFERRED 
EMBODIMENT 

Test SRAM and D-Latch Designs 25 

The test SRAM and latch designs disclosed herein 
were submitted to the MOS Implementation System 
(MOSIS) and fabricated at a 1.6-pm n-well double, 
metal CMOS/bulk foundry. 30 

A schematic diagram of the test SRAM cell is shown 
in FIG. 1. This figure does not include the read/write 
transistors. The pulsed current source is used to model 
an alpha particle strike on drain Dn2 when calculating 
the critical charge of the cell with SPICE. This cell 35 
differs from that of a standard six-transistor SRAM cell 
in three ways: (1) the source of the p-MOSFET, Mp2, is 
connected to an adjustable offset voltage, V 0 , instead of 
Vopto provide a control of the cells critical charge; (2) 
the drain area of n-MOSFET Mn2, Dn2, has been en- 40 
larged by a factor of four over minimum to enhance 
upset rates, thus reducing measurement time; and (3) the 
cell is imbalanced by widening Mn2 over minimum to 
enhance its SEU sensitivity versus V 0 response. The 
dimensions of the MOSFETs and their drains is given in 45 
Table 1. The bloated drain Dn2, is fabricated in a p-type 
substrate doped to 5 X 10 14 cm 3 . The Dn2 junction has 
enhanced charge collection due to funneling compared 
to charge collected by the p-MOSFET drain formed in 
the n-well, which truncates the ion-induced plasma 50 
track. 


TABLE 1 


DEVICE 

L (fim) 

W (/xm) 

AD (jum 2 ) 

Mnl 

1.6 

2.4 

17.92 

Mn2 

1.6 

3.2 

74.88 

Mpl 

3.2 

2.4 

12.16 

Mp2 

3.2 

2.4 

12.16 


In operation all the memory cells are written into a 60 
“sensitive” state where Mn2 is turned ON, connecting 
V 0 to the bloated drain, Dn2. V Q is then lowered from 
Vjm)= 5 V for a period called the stare time. Thereafter 
Vo is returned to Vdd and the cells are read to deter- 
mine the number of upset cells. This cycle is repeated at 65 
different values of V 0 and can be repeated a number of 
times at a given V 0 to improve the resolution of the 
measurement. 


4 

The latch results presented herein were obtained 
from a 64-cell modified transparent D-latch array fabri- 
cated through the same foundry as the test SRAM. In 
FIG. 2 a schematic diagram of this latch is shown with 
the latch in the zero state where reverse biased SEU- 
sensitive drain diodes are represented by said squares. 
The offset voltage, V G was added to sensitize the latch 
to alpha particles. The source of Mp2 is tied to V dd in 
the unmodified latch. The diodes designated with solid 
squares collect the most charge when struck by an ion- 
ized particle because they sit in the p-substrate. The 
layout of the latch cell is shown in FIG. 3. This layout 
includes an inverter that generates enable-bar, and in- 
verters that buffer the Vl-node and the V2-node. The 
sensitive diode areas at the V2-node and V3-node that 
are designated with solid squares in FIG. 2, are outlined 
with a bold line in FIG. 3. These diodes each have an 
area of 69 jim 2 . 

Spice Analysis 

The critical charge, Q c of the test SRAM cell and the 
D-latch cell were determined as a function of V 0 using 
MOSIS supplied parameters and SPICE. The parasitic 
nodal capacitances were modeled by fixed metal and 
polysilicon interconnect capacitances and by the drain 
depletion capacitances given by their areas and periph- 
eries and by the SPICE LEVEL 2 junction capacitance 
parameters. A triangle current pulse with a 1:19 rise:fall 
shape was used to upset these cells. For a given pulse 
height, the transient simulation was evaluated out to 100 
ns where the response was compared to V dd/ 2 to de- 
termine if the cell had flipped. The pulse height was 
adjusted using the binary search algorithm until the 
difference in charge (area under current pulse) between 
successive simulation runs differed by less than 1 fC. 
The resulting critical charge versus V 0 response for the 
V2-node of the test SRAM cell and for the average of 
the V2-node and V3-node for the D-latch is shown in 
FIG. 4. These results were found to be invariant with 
current pulse widths up to 500 ps. Since the alpha parti- 
cle pulse width is about 200 ps, the response of these 
circuits exceeds that of the alpha particle pulse width. 

As V 0 decreases the critical charge of the SRAM cell 
and D-latch decreases approximately linearly and goes 
to zero at V 0 = V 05 . For the SRAM cell V 0 sS~ L8 V and 
for the latch Y 05 l=2A V. DC SPICE simulations of 
these bistable circuits show that the sensitive state and 
metastable state approach each other as V 0 decreases 
until they are equal at V 0 =V 0S . For V 0 ^ V 05 the circuit 
spontaneously goes to the flipped state. The circuit at 
this point is no longer bistable. V 05 is a function of the 
threshold voltages and geometries of the n- and p-MOS- 
FETs. Process induced dispersion in the MOSFET 
parameters gives rise to a dispersion in V os for an array 
of cells. The slopes of the critical charge curves in FIG. 
4 are defined by their upset capacitances, Cc/=8Q C /- 
8V 0 . The upset capacitance of the SRAM cell is 
Ct/s=56 fC/V and of the latch cell is C^=850 fC/V. 
SPICE is a well-known computer program for simulat- 
ing and analyzing circuits which, for example, is de- 
scribed in detail in the text entitled “SPICE: A Guide 
To Circuit Simulation And Analysis Using PSPICE”, 
second edition by Paul W. Twimenga, copyright 1992, 
Prentice Hall. 

Test SRAM Alpha-Particle Response 

The response of the SRAM to alpha particles from a 
4.6 jLtCi Am-241 source placed 0.508 cm above the 
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SRAM in air is shown in FIG. 5. Also shown, is the 
spontaneous response which was measured without the 
alpha-particle source present. As V 0 is lowered from V 0 , 
the SRAM begins to upset at V 0 <3.2 V in a region 
termed the tail region. As V 0 is further lowered, the 5 
upset rate increases more slowly; this region is termed 
the periphery hit region. Finally, as V 0 reaches V os = 1.8 
V, the upset rate increases very rapidly in what is 
termed the spontaneous flip region. 

In the periphery hit region, cells upset when the alpha 10 
particle hit the Dn 2 drain or come close to, but miss, the 
Dn2 drain. This effect can be explained by a delayed 
fleld-funneling effect and to diffused charge that is col- 
lected by Dn2. Since the cell critical charge decreases 
as V 0 decreases, the alpha particles can strike further I 5 
from the Dn2 depletion edge and still upset the cell. 

The upset rate R, displayed in FIG. 5, was obtained 
by dividing the average measured number of upsets 
that occurred in a stare time t by t, or R=N a /t, where 


1 k 1 ( k y 

25 

and the number of samples, k— 200. The stare time in 
this case was fixed at 2 seconds. The spontaneous upset 
data was also divided by t so that both responses satu- 
rate at N//t=2048, where N* is the number of ceils in 
the memory cell (4096 in this case). The response of the 3 
SRAM cells is given by the detector equation: 


6 

The tail region response shown in FIG. 5 was 
smoothed and differentiated using a twenty-five point 
least squares difference algorithm. This differential re- 
sponse, shown in FIG. 7, is composed of a number of 
peaks plus the onset of the peripheral hit response. 
These peaks are due to the selective degradation of the 
alpha-particle energies as the particles pass through the 
various regions found in the bloated drain Dn2. The 
major peak, shown in FIG. 7, is due to the alpha parti- 
cles that pass through the thickest regions; those cov- 
ered with metal 2 (M2). This M2 region is 38.2 pm 2 . In 
the region covered with M2, the alpha particles lose 
more energy than in other regions so that they have the 
highest linear energy transfer (LET). 

A schematic view of the path of the alpha particle 
through the silicon is shown in FIG. 8. The thickness 
X4— X2 is obtained by first converting the full- width- 
half-maximum of the peak in the response curve shown 
in FIG. 7, V 0 f, to the voltage straggle, V 0(r , of the alpha 
particles at the end of the collection region in the silicon 
using V 0(r —Vof/ 2.35. Then the voltage straggle is con- 
verted to energy straggle using 

E (r =(C u /K)V 0<r (5) 

Finally, X4—X2— 13.58 pm was calculated using the 
range-straggle tables for alpha particles in silicon. 

The silicon collection layer thickness 8X4 was deter- 
mined from the offset voltage shift 8V 0 =V op — 
V 0 s— 0.98 V. This was converted to an energy shift 
using 


- or <t>(N t - N) 


( 2 ) 


whose homogeneous solution with initial condition at 
t=0 of no tripped cells (N=0) is: 

AT^iVKl-c- 0 ^) (3) 

where cr is the device cross section (Dn2 drain area), 
and O is the total integrated flux whose particles deposit 
charge exceeding the critical charge Q c required to trip 
a cell. 

The upset capacitance Cu, the spontaneous upset 
voltage V 0 s, and the hole-electron pair charge-energy 
factor for silicon K=44.2 fC/MeV enable one to cali- 
brate the offset voltage V 0 axis in terms of the critical 
energy E c , by 

E c ^(Cu/K)(V 0 -V 0S ), (4) 


hE^(C u /K)V 0 (6) 

35 8X4=6.09 pm was calculated using the energy range 
tables for alpha particles in silicon. Finally the over- 
layer thickness, 8X3 was calculated by 

6X3 =(X4-X 2 )- 8X4 =7.59 pm. 

4Q D-Latch Alpha Particle Response 

The alpha-particle response of the latch array was 

obtained in the same manner as for the SRAM except 
that the Am-241 source was placed 0.558 cm above the 
latch chip in air. The results, shown in FIG. 9, include 
45 the spontaneous data obtained before and after exposing 
the part to alpha particles. The shift observed in these 
responses is due to total dose effects such as the charg- 
ing of gate oxides. The average of the pre- and post- 
radiation spontaneous data was used in the analysis of 
50 this part. 

The dispersion in the latch alpha-particle response, 


shown in FIG. 9, is primarily that of the spontaneous 
or for this SRAM, E C =1.267(V 0 — 1.8) MeV. response because the offset voltage dispersion due to 

FIG. 6 illustrates the relationship between upset rate the alpha-particle energy spectrum is much smaller than 

R, the SRAM spontaneous density n, the source spectra 55 the dispersion in the spontaneous response. This is a 
3>, and the SPICE derived critical charge curve. The consequence of the reduced energy detection sensitivity 


spontaneous density, n=dR/dV 0 , is the derivative of 
the spontaneous response data shown in FIG. 5, and 
was found to be Gaussian with a mean and standard 
deviation of V^— 1.8 ±0,011 V. The area under the 
n(V 0 ) curve is equal to N ? . The exponential periphery 
hit response is not shown in FIG. 6 because it was sub- 
tracted from the measured R. In the analysis to follow, 
the spontaneous density, n, is approximated as an im- 
pulse function. This is justified for the SRAM data 
shown in FIG. 5 because the slope of n/t=4096 cells/34 
mV sec is much steeper than the SRAM alpha-particle 
response. 


of the latch (K/Cul~ 0.052 V/MeV for the latch ver- 
sus K/C(/s=0.789 V/MeV for the SRAM). The disper- 
sion in the latch alpha-particle response, Vos =34 mV, 
was found by fitting the data to a Gaussian function. As 
shown in FIG. 9 the shift in the offset voltage, SV 0 , was 
obtained between the points on the average spontane- 
ous curve and the alpha-particle curve where, on the 
average, half of the latch cells are flipping (N ? /2=32). 
Using the SPICE derived upset capacitance for the 
latch, C(7£=850 fC/V, the collected charge required to 
upset a latch cell is 8Q C £=50.2 ±28.9 fC which agrees 
with that obtained with the SRAM, Q c s- The dispersion 
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in QcLj however, is large compared to the dispersion in 
QcS . Since Q c l— QcS, it is valid to apply the collection 
depth dX4 and the over layer thickness obtained with 
the SRAM to the analysis of the latch. 

For the latch, the LET threshold for Vdd= 5 V oper- 5 
ation can be predicted from the SPICE derived upset 
capacitance, Cul , the measured mean spontaneous off- 
set voltage, Vos, and the collection depth obtained from 
the SRAM alpha-particle response, 8X4: 

10 

, „ C UL (V DD • V os ) „ , MeV cm 2 ( 7 ) 

LET — K^m — = 34 ' 

where, Cra,=850 fC/V, V os =2.49 V, K=44.2 
fC/MeV, p=2320 mg/cm 3 for Si, and I 5 
8X4=6.09 X 10~ 4 cm. 

D-Latch Heavy Ion Response 

Heavy ion tests were performed at the Brookhaven 
National Laboratory on the latch array with 20 
Vx>d=V 0 = 5 V. These measurements were made by 
gating the ion beam on and then off when a fluence, 
F=1.13X10 7 ions/cm 2 cos 0 was reached, where 0 is 
the angle of incidence. The stare time was set so that an 
average of two upsets would occur each test cycle. This 25 
assured that essentially all N,= 64 latch cells in the array 
were susceptible to upset and allowed the effective 
cross section per cell to be calculated by: 

N (8) 30 

" N t F cosd 

where N is the total number of upsets that occur while 
the ion beam is gated on. 

Using the average over layer thickness, 8X3 =4.3 pm, 35 
and the collection depth, 8X4= 6.1 pm obtained from 
the SRAM alpha particle measurements, the energies of 
the ions used in the test were computed from the TRIM 
at the entry point of the Si collection layer, X3, and at 
the exit point of the Si collection layer, X4. The average 40 
LET of the ion in the collection layer was then calcu- 
lated by: LET=[E(X3)— E(X4)]p8X4. These results are 
shown in Table 2 and in FIG. 9. 


TABLE 2 


ION 

PLOT 

SYMBOL 

ANGLE, 6 
(degrees) 

cr 

Oxm 2 ) 

LET 

(MeV/cm 2 mg) 

1-127 

□ 

0 

232 

58.8 

Br-79 

■ 

45 

209 

57.3 

Br-79 

O 

35 

144 

49.5 

Br-79 

• 

0 

63.8 

40.3 

Fe-56 

X 

45 

33.0 

37.5 


A test SRAM, designed to be upset-sensitive (LET 
soft) and having an electrically adjustable critical 
charge has been characterized using alpha particles 55 
using an inexpensive laboratory setup. The SRAM 
alpha particle results and SPICE simulation results fa- 
cilitated in predicting the LET threshold of a standard- 
cell D-latch. The predicted LET =34 MeV cm 2 /mg 
compares favorably with heavy ion test value of 41 60 
MeV cm 2 /mg. 

Preferred Analysis Method 

A preferred method for calibration of the SRAM 
detector will now be described in conjunction with 65 
FIGS. 11 to 16. The procedure requires determining the 
detector’s overlayer thickness, 8X3, and collection 
depth, 8X4. These were determined using 0.55 and 1.0 
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MeV protons from the Caltech Tandem Van de Graaff. 
The analysis requires knowledge of the proton charge 
deposition versus range and this was calculated using 
particle range physics from TRIM. The validity of the 
analysis was verified with energy straggling measure- 
ments. Finally, the SPICE circuit simulation program 
was used to compute the relationship between the 
charge deposited in the memory cell and the cell offset 
voltage. 

The critical charge, Q c , of the SRAM cell was deter- 
mined as a function of V 0 using MOSIS supplied SPICE 
parameters. The parasitic nodal capacitances were mod- 
eled by fixed metal and polysilicon interconnect capaci- 
tances and by the drain depletion capacitances using 
their areas and peripheries. The SPICE simulation used 
the level-2-model. A triangle current pulse with a 1:19 
rise:fall shape was used to upset these cells. The location 
of the pulse generator in the memory cell is shown in 
FIG. 1. For a given pulse height, the transient simula- 
tion was examined at 100 ns where the response was 
compared to V dd / 2 to determine if the cell had flipped. 
The current pulse height was adjusted using a binary 
search algorithm until the difference in charge (area 
under current pulse) between successive simulation runs 
differed by less than 1 fC. These results were found to 
be invariant with current pulse widths up to 500 ps. 
Since the proton width is about 200 ps, the response of 
these circuits exceeds that of the proton current pulse. 

SPICE simulations provide the Q c versus V 0 curve 
shown in FIG. 11 which is well approximated by the 
straight line relationship: 

Q c =Cu(V o -V 0SH ) 

where V 0:fft =1.8 V is the mean offset voltage in the 
spontaneous flip range and the slope is the upset capaci- 
tance, C u = dQ c /d V 0 =56 fC/V. 

In the SRAM test sequence, all memory cells are 
written into the sensitive or stare state where Mn2 is 
turned OFF and Mp2 is turned ON which connects 
to the bloated drain, Dn2; see FIG. 1. V 0 is then low- 
ered from 5 V for a period called the stare time which 
can last a few seconds to several days. Thereafter, V 0 is 
returned to 5 V and the cells read to determine the 
number of upsets. This cycle is repeated for different 
values of V 0 . 

Test results are shown in FIG. 12 for protons and 
alpha particles and compared against the spontaneous 
flip response. The overall behavior indicates that at 
high V 0 the response curve is determined by particle 
energy straggling. At intermediate V 0 , the response is 
determined by the collection of particles outside the 
area of drain Dn2. This is called the peripheral hit re- 
gion of the curve. Finally, at low V c , the cells flip spon- 
taneously. 

Layer Thickness Analysis 

The following analysis is used to determine the 
SRAM detector Dn2 overlayer material thickness, 8X3, 
and charge collection region thickness, 8X4. A sche- 
matic view of these layers is shown in FIG. 13. 

The analysis requires the determination of offset volt- 
age, V 0 p, at the peak of the particle upset distribution. In 
this technique the particle upset distribution is Gaus- 
sianly distributed as seen in the cumulative probability 
plot shown in FIG. 14. The analysis follows from the 
SRAM detector equation: 
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dN/dt= cr$(N t — N) 

where cr is the area of Dn2,<f> is the flux, N is the mea- 
sured number of flipped cells, and Nr is the total number 5 
of cells in the SRAM. Evaluating for N->0 

dN/dt N-*0=R o =cr<\>Nt 

The particle flux at diode, Dn2, is described by the 10 
Gaussian or normal distribution expressed by the error 
function: 

4> = 4>mU ~ emy op - V op p)/V op <r V7]}/2 15 

where 4> m is the maximum particle flux entering the 
SRAM, V 0 p is the offset voltage in the vacinity of the 
upset peak, V opfl is the mean or peak value, and V op(r is 
the standard deviation for the upset distribution. Note 20 
that (f> = (p m /2 at V 0/ ,==V W which states that only half 
the particle flux entering the SRAM can upset cells. 

Combining the above equations leads to the following 
expression for the number of flipped cells in the vacinity 
of the particle upset distribution: 25 

N = Ro/<r4> m = N,{1 - eWVop - V op n)/V op o-M7]}/2 

30 

The probability distribution given in percent and 
plotted in FIG. 14 was determined from: 

P{V oi > V op )= KXHiV— 0.5)/iV, 

where N is the number of flipped cells at V^and for this 35 
SRAM N,=4096 cells. 

The calculation of N can take two approaches. The 
first approach is used when the experimenter knows the 
beam flux, Q m . In this case, N^R^/crcf)™ where or =68.8 
jam for Dn 2. The second approach is used when the 
experimenter knows the beam fluence, F m , and the 
number of upsets, N/, observed during the time the 
beam is on. That is N=N/crF m . This approach is useful 
when the flux is variable and the fluence can be moni- .. 
tored. 45 

The data points in FIG. 14 allow an analysis of the 
Gaussian nature of the energy dispersion of the particles 
as they lose energy in the silicon. At the peak of the 
upset distribution, at the 50 percent point, only half of 5Q 
the particles can deposit sufficient charge to flip the 
cells. This is defined as Q cp and is given by: 

Qcp=C u (V opp - Vosp) 

where Q cp values for the four curves shown in FIG. 14 55 
are listed in Table 3. 

The SRAM overlayer thickness, 8X3, was deter- 
mined using a 0.55 MeV proton beam which stopped 
within the collection layer. This is crucial for it means 
that all the charge deposited in the collection layer, 
8X4, is collected by the diode Dn2. The charge depos- 
ited by the 0.55 MeV proton beam is shown in FIG. 15. 
This curve was plotted using TRIM. The analysis for 
8X4 begins by determining V opp for the 0.55 MeV pro- 65 
ton upset data shown in FIG. 14; the results are listed in 
Table 3. Then Q c/7 is determined from FIG. 11. Finally, 
the Q cp value is subtracted from the End of Range, 
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EOR, of the charge deposition curve shown in FIG. 15. 
Thus 8X3=4.32 jam. 

The charge collection thickness, 8X4, was measured 
with 1.0 MeV proton beam which had a range greater 
than 8X3 +8X4. The charge deposited by the 1.0 MeV 
proton beam is shown in FIG. 16. The analysis begins 
by determining V # for the 1.0 MeV proton upset data 
shown in FIG. 14. Then Q cp was determined from FIG. 
11 . Finally 8X4 was determined by adding Q cp to the 
overlayer charge as determined from 8X3; the tech- 
nique is illustrated in FIG. 16. Thus 8X4=6.64 jam. 

The charge collection depth, 8X4, was also deter- 
mined for a 4.7 MeV alpha particle and a collection 
depth of 6.33 jam was determined. This value is slightly 
smaller than the collection layer thickness determined 
for the 1.0 MeV proton. This is due to the fact that 
heavy particles have a shorter range. 

The linear energy transfer, LET, in MeV -cm 2 /mg for 
the particles can now be calculated using: 

LET = C u { V opp — V 0S}l )/KpdX, 4 

where for silicon K=44.2 fC/MeV, p= 2320 mg/cm 3 , 
C u is in fC/V, V 0 is in Volts, and 8X4 is in cm. LET 
values are listed in Table 3. The LET value for V w =5 
V is 2.88 MeV cm 2 /mg for an estimated 8X4=6.00 pm. 
The conversion factor K is determined from: 

K= 1.6 X 10 “ 19 (CouI./electron)X 10 15 
(/Coul./Coul.)X 10 6 (eV/MeV)/3.62(eV/eh-pm). 

Layer Thickness Dispersion Analysis 

The result of calculating the errors in 8X3 and 8X4, 
that is, 8X0- and SXo-, is shown in Table 3 . It shows that 
as the 0.55 and 1.0 MeV protons pass through the vari- 
ous regions that the energy dispersion increases. This 
leads to an uncertainty in the thickness of each layer. 
The thickness dispersion is calculated from the energy 
dispersion. For the over-layer the depth dispersion is: 

SX 3<r = SE 3 X <T -K-ddX 3 /dQ d 

where 8E3X<ris the energy dispersion due to over-layer 
thickness variations and from FIG. 15, dQj/d8X3=3.85 
fC/jam. The collection layer depth dispersion is: 

8X 4o .=dE 4 X <r -K-ddX 4 /dQ d 

where 8E4XO* is the energy dispersion due to collection 
depth variations and from FIG. 16, dQ^/d8X4=3.125 
fC/jam. 

The evaluation of 8X30- requires an evaluation of 
8E30-. Since the 0.55 MeV protons stop in the charge 
collection region, the energy dispersion is determined 
by the Au scattering foil and the Si over-layer and not 
the Si collection layer. The energy dispersion for the 
0.55 MeV protons is evaluated at the peak of the proba- 
bility distribution as defined at the 50 percent point in 
FIG. 14, it is calculated from E 0/JO - (0.55) = (C u /K.yV op<r 
(0.55). This energy dispersion consists of the following 
components: 

E op<T (0.55) = [(E 0S{r ) 2 + (£ /or ) 2 + (S£ 3 o-)¥ 

where E 05(T = (C w /K) V 05Cr is the instrument function 
energy dispersion, E/o- is the source Au scattering foil 
energy dispersion, and 8E3o- is 
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where SEsSo- is over-layer energy dispersion due to 
particle straggling and 8E3Xcr is the over-layer energy 
dispersion due to over-layer thickness fluctuations. The 5 
E/cr=0.012 MeV was determined from Rutherford scat- 
tering theory for a 1.2-/xm thick Au foil. 8E3So-=0.009 
MeV was determined from Rutherford scattering the- 
ory for a 4. 32- jam thick over layer. SEsXo-is now calcu- 
lated from the above two equations and 8X30-— 0.19 jam 10 
was then determined. 

The evaluation of 8X40- requires an evaluation of 
8 E 4cr . Since the 1.0 MeV protons pass through the 
charge collection region, the energy dispersion is deter- 
mined by the Au scattering foil, over-layer and collec- 15 
tion layer. The energy dispersion for the 1.0 MeV pro- 
tons is evaluated at the peak of the probability distribu- 
tion as defined at the 50 percent point in FIG. 14 ; it is 
calculated from E 0 pcr (1.0)=(C a /K) V 0 p<r (1*0). This 
dispersion consists of the following components: 20 

Eopo- (1*0)= [(^<r) 2 +(£/ ( r) 2 +(S£3cr) 2 +(6^) 2 ] i 


where 8E4cr is 

SE 4 <r= [(SEaSo-) 2 HSEaXv) 2 ]* 


25 
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1. A method of predicting the SEU susceptibility of a 
SRAM; the method comprising the steps of: 

a) providing an enhanced SEU sensitivity test SRAM 
having an electrically adjustable critical charge; 

b) determining the critical charge as a function of 
offset voltage to ascertain the upset capacitance 
C us of said test SRAM; 

c) exposing the test SRAM to a first known level of 
alpha particle energy to ascertain the overlayer 
thickness of said SRAM, said first known level 
being selected to stop the alpha particle energy in 
the collection layer of said SRAM; 

d) calculating the overlayer thickness 8X3 of said test 
SRAM based upon the charge deposited by said 
first known level; 

e) exposing the test SRAM to a second known level 
of alpha particle energy to ascertain the charge 
collection thickness of said SRAM, said second 
known level being selected so alpha particle energy 
stops beyond the collection layer; 

f) calculating the collection layer thickness 8X4 of 
said test SRAM based upon charge deposited by 
said second known level and the calculated thick- 
ness of said overlayer; and 

g) calculating the linear energy transfer (LET) for 
said test SRAM: 


where 8E4S0- is collection layer energy dispersion due 
to particle straggling and 8E4Xo- is the energy disper- 
sion due to collection layer thickness fluctuations. The 30 
8E4 So-= 0.011 MeV was determined from Rutherford 
scattering theory for a 6.64-ju,m thick collection layer. 
8E4 Xo-= 0.025 MeV is now calculated from the above 
two equations and 8X40- =0.3 6 /xm was then deter- 
mined. 

TABLE 3 


8X3, 8X4, and LET Values 


PARTICLE 

Qcp 

fC 

V opp 

V 

V opcr 

V 

8X3 

pm 

8X4 

pm 

LET 

MeV* 

cm 2 /mg 

0.55 MeV 
PROTON 

10.47 

1.987 

0.021 

4.32 

- 

— 

1.00 MeV 
PROTON 

16.18 

2.089 

0.030 

4.32 

6.64 

0.24 

4.70 MeV 
ALPHA 

50.40 

2.66 

0.143 

4.32 

6.33 

0.77 

SPONTANEOUS 

178.00 

5.00 

v 0 ^ 

1.80 

, §2 

I 

> 0 

4.32 

6.00 

2.88 


Having thus described an exemplary embodiment of 
the invention, what is claimed: 


T FT _ C US (V 0 p\i - V os p) 

“ Kp8X 4 

where 

V 0 pp~ the peak value of offset voltage at maximum 
particle upset distribution; 

V 05/A =the mean offset voltage in the spontaneous flip 
range; 

K=the hole-election pair charge-energy factor for 
silicon; 

p=the density of silicon; 

Cus — upset capacitance for SRAM. 

40 2. The method recited in claim 1 wherein in step a) 

said SRAM is made more sensitive to incident-ion- 
induced SEU by providing at least one pull-down field 
effect transistor having a bloated drain surface area and 
at least one pull-up field effect transistor having a 
45 source connected to an offset voltage. 

3. The method recited in claim 1 wherein said offset 
voltage is determined by measuring the offset voltage at 
which about one-half of the total number of latch cells 
are upset spontaneously without deliberate ion expo- 
50 sure. 
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